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If} . Abstract 

By using thermal QCD sum rules to investigate the 8 induced neutron 
, electric dipole moment d n , we have examined the behaviour of broken CP 

symmetry at finite temperature. We find that, below the critical tempera- 
ture, the ratio | ~ | slightly decreases but survives at temperature effects, 
Ph, implying T nonrestoration of CP-invariance Q . 



INTRODUCTION 



To establish connection between particle physics and cosmology, it is essential 
to study the behaviour of symmetry breaking in the early universe, i.e. at high 
temperature. The CP symmetry is, without doubt, one of the most fundamental 
symmetries in nature. It is intimately related to theories of interactions between 
elementary particles and represents a cornerstone in constructing grand unified and 
super symmetric models. However, its breaking still carries a cloud of mystery in 
particle physics and cosmology since it is necessary to explain baryogenesis and is 
required by theories with domain walls. 

According to the CPT theorem, CP violation implies T violation. The latter is 
tested through the measurement of the neutron electric dipole moment (NEDM) 
d n . The upper experimental limit gives confidence that the NEDM can be another 
manifestation of CP breaking. In the Standard Model, CP violation is originated 
from two sources: the first source, which appears in the electroweak sector, is 
parametrized by a single phase in the Cabbibo-Kobayashi-Maskawa (CKM) quark 
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mixing matrix [[|. The other source is due to the the so called ^-term of QCD. In 
fact, the the QCD effective lagrangian contains an additional CP-odd four dimen- 
sional operator embedded in the following topological term: 



U = (1) 

8"7T 

where is the gluonic field strength, G^ u is its dual and a s is the strong 
coupling constant. The G fJiU G' J ' u quantity is a total derivative, consequently it 
contribute to the physical observables only through non perturbative effects. The 
NEDM is related to the #-angle by the following relation: 



e m q - 2.7 x 10- 16 ^ § 
4 W n { Wj ~ { 5-2x 10^9 § (2) 

and consequently, according to the experimental measurements d n < 1.1 x 
10~ 25 ecm ||, the 9 parameter must be less than 2 x 10 -10 ||. The well known 
strong CP problem consists in explaining the smallness of 9. In this regard, several 
scenarios were suggested. The most popular one is the Peccei and Quinn 0, in 
which 9] is identified to a very light pseudo goldstone boson called the axion. This 
particle arises from the spontaneous breakdown of a global Ua(1) symmetry and 
may well be important to explain the puzzle of dark matter providing a peace of 
information on the missing mass of the universe ||. 

Our aim in this work is to study the thermal behaviour of the CP symmetry 
breaking and the temperature effects on its restoration. This is motivated by the 
possibility to restore some broken symmetries by increasing the temperature. 

This paper is organized as follows: Section 2 is devoted to the calculations of 
the NEDM induced by the 9 using QCD sum rules. In section 3, we show how one 
introduces temperature in QCD sum rules calculations. The last section is devoted 
to a discussion and qualitative analysis of the thermal effects on the CP symmetry. 



NEDM from QCD sum rules 

Since the NEDM has essaentially a non perturbative nature, one should any how 
take into account some effects which escape the perturbative treatement. One 
such approach, dealing with the strong strong coupling regime and based on first 
principales of the theory is QCD sum rules. Such approach has been applied 
successfully to the investigation of hadronic properties at low energies, particularly 
to certain baryonic magnetic form factor. In order to derive the NEDM through 
through this approach || |TD|, we consider a lagrangian containing the following P 
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and CP violating operators: 

L p ,cp = -e q m. J] q f i l5 q f + B-^-G^GT. (3) 

9 q and # are respectively two angles coming from the chiral and the topological 
terms and m* is the quark reduced mass given by m* = ™ u +^ d ■ Then, as usual, 
we start from the two points correlation function in QCD background with a non- 
vanishing 9 and in the presence of a constant external electomagnetic field F^ v 

Tl(q 2 ) = i / d 4 xe iqx < 0\T{r](x)f](0)}\0 > 9)F . (4) 
where rj(x) is the neutron interpolating current fl~2|| : 



rj = 2e abc {(d a C'j 5 u b )d c + f3(d a Cu b )^d c } 1 (5) 

and (3 is a mixing parameter. To select the appropriate Lorentz stucture, H(q 2 ) is 
expanded in terms of the electromagnetic charge as: 

II(g 2 ) = n (0) (g 2 ) + en (1) (g 2 ,FH + O(e 2 ). (6) 

The first term rT°)(g 2 ) is the nucleon propagator which includes only the CP-even 
parameters, while the second term H^\q 2 , F^ u ) is the polarization tensor which 
may be expanded through Wilson OPE as: C n < 0|gTg|0 >e )F , where Y is an 
arbitrary Lorentz structure and C n are the Wilson coefficient functions calculable 



in perturbation theory [14, 13 1. From this expansion, we keep only the CP-odd 



contribution part. The electromagnetic dependence of these matrix elements is 



determined in terms of the magnetic susceptibilities k, x an d £, defined as [ |13|| : 



< 0|gV^g|0 > F = xe g F^ < 0\qq\0 > (7) 
g < 0\qG^q\Q > F = Ke q F^ < 0\qq\0 > (8) 
2g < OlqG^qlO > F = fe^F"" < 0|gg|0 > (9) 

Besides, by considering the anomalous axial current, one obtains the following 
9 dependence of < 0|gTg|0 >e matrix elements 0: 

m q < 0|gTg|0 > e = imj < 0\qTq\0 > +0{m 2 q ) (10) 

where m q and m* are respectively the quark and reduced masses. 0(m 2 ) con- 
nection is negligible since m v >> m w . 

Putting altogether the above ingredients and after a straightforward calculation 
10 1, the following expression of U^(q 2 , F^ u ) for the neutron is derived: 
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n(-g 2 ) = -^<0|^|0>{F ( T,g}[x( / 9 + l) 2 (4e d -e u )ln(^ 
-4(0 - l) 2 e„(l + ^(2 K + 0)(ln(^) - 1) ' 



2 



4 tf R ~q 2 

^((4/5 2 - 4/5 + 2)e d + (3f3 2 + 2/5 + l)e u )^...], (11) 
2 — 



where # = 9 + # g is the physical phase and q = <3^7 M . 
The QCD expression (2.7) is confronted to the phenomenological parametrisation 
U Phen (—q 2 ) written in terms of the Neutron hadronic properties. The latter is 
given by: 

a ^ ) ^ [ p aMj ^ + w ±_ + ... ) , (12) 

where m n is the neutron mass, e q is the quark charge. A and A 2 , which originate 
from the phenomenological side of the sum rule, represent respectively a constant 
of dimension 2 and the neutron coupling constant to the interpolating current r](x). 
This coupling is defined via a spinor v as < 0|?7(x)|n >= Xve 01 " 75 . 



Thermal NEDM sum rules 

The introduction of finite temperature effects may provide more precision to the 
phenomenological values of hadronic observables. Within the framework of QCD 
sum rules, the Temperature evolution of the correlation functions manifests itself in 



the thermal average of the Wilson operator expansion |15|, [LSj] . Hence, at relatively 
low temperature, the system can be regarded as a non interacting gas of bosons. 
In this approximatyion, the Thermal dependance of the vacuum condensates can 
be written as : 

< & > T =< & > +/ < 7r(p)|OV(p) > n B (±) (13) 

where e = Jp 2 + m 2 , Ub = the Bose- Einstein distribution and < O l > is the 

standard vacuum condensate (i.e. at T=0). In this approximation, we only kept 
the pion contributions, since in the low temperature region, the effects of heavier 

-mp 

resonances (r = K, rj, ..etc) are dumped by their distibution functions ~ e~r~[|17|. 
To compute the pion matrix elements, we apply the soft pion theorem given by: 

1 77? 2 

< ApWHp) >= —r 2 < 0\[QI [Q a 5 ,OW > +0(-f ), (14) 

J 7T 
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where A is a hadron scale and Q5 is the isovector axial charge defined by: 

Qt = J (Pxq(x)^ 5 —q(x). (15) 
Direct application of the above formula to the quark and gluon condensates shows 



the following features [16|, [17] : 



^i) Only < qq > is sensitive to temperature. Its behaviour at finite T is given by: 

< qq > T ~ (1 - < m >, (16) 



where <p(T) = ^B(^) with B(z) = 4, f™ dy^^f- and f n is the pion decay 
constant (f n ~ 93MeV). The variation with temperature of the quark condensate 
< qq >t results in two different asymptotic behaviours, namely: 

T 2 m 
< qq > T ~ (1 - — ) < qq > for < 1 (17) 

lnm n T 2 -rnn _ 
< qq > T ~ (1 - ^ ^g^e t ) < qq > for — > 1 (18) 

(ii) The gluon condensate is nearly constant at low temperature and a T de- 
pendence occurs only at order T 8 . 

The determination of the ratio ~ sum rules at non zero temperature is now 
easily performed by applying Borel operator to both parametrisation of the Neutron 
correlation function shown in Eqs. (2.7) and (2.8). Then finite temperature effects 
are introduced via the procedure discussed above. Finally, by invoking the quark- 
hadron duality, we deduce the final sum rules of the 9 induced NEDM at finite 
temperature: 



cL , M 2 m* 1 . <p(T) s r , £ , Ml 

T (T) = - UP WWJT) { 1 8 } < ™ > te(4e»-e i )-^(4e„+8e J )]^ 

(19) 

where M represents the Borel parameter. The single pole contribution entering the 
sum rules via the constant A has been neglected, as suggested in 0. 
The value of (3 has been set to 1 in (3.5). This choice is more appropriate for 
us since it suppresses the infrared divergences. The T- evolution of the coupling 
constant and the mass of the neutron were determined from the thermal nucleon 



sum rules [|T7 |. 



Within the dilute pion gas approximation, Eletsky has shown that the contri- 



bution induced the pion- nucleon scattering has to be considered ||19||. It enters the 



nucleon sum rules through the coupling constant g n NN, whose values lie within the 
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range 13.5-14.3 @ . 

Numerical analysis is performed with the following input parameters: the Borel 
mass has been chosen within the values M 2 = 0.55 — 0.7GeV 2 which correspond to 
the optimal range (Borel window) in the % sum rule at T = ]TU|] . For the x an d 
£ susceptibilities we take x = -5-7 ± Q.QGeV' 2 [pTJ and f = -0.74 ±0.2 p|. As 



to the vacuum condensates appearing in (3.5), we use their standard values |TT 



Analysis and Conclusion 

We have established the relation between the NEDM and 9 angle at non zero 
temperature from QCD sum rules. We find that the behaviour of the ratio^ 1 is 
connected to the thermal evolution of the pion parameters /„-, m n and of g n NN- 
By analysing the ratio as a function of T in the region of validity of thermal sum- 
rules [0,T C ], we learn that | | decreases smoothly with T (about 16% variation 
for temperature values up to 200 MeV) but survives at finite temperature. This 
means that either the NEDM value decreases or 9 increases. Consequently, for a 
fixed value of 9 the NEDM decreases but it does not exhibit any critical behaviour. 
Furthermore, if we start from a non vanishing 9 value at T = 0, it is not possible 
to remove it at finite temperature. We also note that | % | grows as M 2 or x 
susceptibility increases. It also grows with quark condensate rising. However this 
ratio is insensitive to both the £ susceptibility and the coupling constant g n NN- We 
notice that for higher temperatures, the analysis of | ~ \= f(^r) exhibits a brutal 
increase justified by the fact that for temperatures beyond the critical value T c , 
at which the chiral symmetry is restored, the constants and g n NN become zero 
and consequently from Eq(3.5) the ratio % behaves as a non vanishing constant. 
The large difference between the values of the ratio for T < T c and T > T c may 
be a consequence of the fact that other contributions to the the spectral function 
have ben neglected, like the scattering process N + ir — > A. These contributions 
which are of the order T 4 , are negligible in the low temperture region but become 
substantial for T > T c . Moreover, this difference may also originate from the use 
of soft pion approximation which is valid essentially for low T (T < T c ). Therefore 
it is clear from this qualitative analysis, which is based on the soft pion approx- 
imation, that temperature does not play a fundamental role in the suppression 
of the undesired (9-term and hence the broken CP symmetry is not restored fl24 |. 



Indeed, some exact symmetries can be broken by increasing temperature 24] . 
The symmetry non restoration phenomenon, which means that a broken symmetry 
at T=0 remains broken even at high temperature, is essential for discrete symme- 
tries, CP symmetry in particular. Indeed, the symmetry non restoration allows us 



to avoid wall domains inherited after the phase transition p5[ and to explain the 



baryogenesis phenomenon in cosmology pfifl . Furthermore, it can be very useful 



for solving the monopole problem in grand unified theories p7 . 
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